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Abstract 

We propose a novel approach to determine the leading hadronic corrections to the muon g-2. It consists 
in a measurement of the effective electromagnetic coupling in the space-like region extracted from Bhabha 
scattering data. We argue that this new method may become feasible at flavor factories, resulting in an 
alternative determination potentially competitive with the accuracy of the present results obtained with the 
dispersive approach via time-like data. 


1. Introduction 

The long-standing discrepancy between experi¬ 
ment and the Standard Model (SM) prediction of 
a^, the muon anomalous magnetic moment, has 
kept the hadronic corrections under close scrutiny 
for several years m- In fact, the hadronic uncer¬ 
tainty dominates that of the SM value and is com¬ 
parable with the experimental one. When the new 
results from the g-2 experiments at Fermilab and 
J-PARC will reach the unprecedented precision of 
0.14 parts per million (or better) [SHZ]) the uncer¬ 
tainty of the hadronic corrections will become the 
main limitation of this formidable test of the SM. 

An intense research program is under way to 
improve the evaluation of the leading order (LO) 
hadronic contribution to due to the hadronic 
vacuum polarization correction to the one-loop di¬ 
agram ElE], as well as the next-to-leading order 
(NLO) hadronic one. The latter is further divided 
into the 0{a^) contribution of diagrams contain¬ 
ing hadronic vacuum polarization insertions nni, 
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and the leading hadronic light-by-light term, also 
of 0{a^) [5J[TT1[T2]. Very recently, even the next- 
to-next-to leading order (NNLO) hadronic contri¬ 
butions have been studied: insertions of hadronic 
vacuum polarizations were computed in m, while 
hadronic light-by-light corrections have been esti¬ 
mated in El- 

The evaluation of the hadronic LO contribution 
involves long-distance QCD for which pertur¬ 
bation theory cannot be employed. However, us¬ 
ing analyticity and unitarity, it was shown long 
ago that this term can be computed via a disper¬ 
sion integral using the cross section for low-energy 
hadronic e“'"e“ annihilation m- At low energy this 
cross-section is highly fluctuating due to resonances 
and particle production threshold effects. 

As we will show in this paper, an alternative 
determination of can be obtained measuring 
the effective electromagnetic coupling in the space¬ 
like region extracted from Bhabha (e+e” —)■ e+e”) 
scattering data. A method to determine the run¬ 
ning of the electromagnetic coupling in small-angle 
Bhabha scattering was proposed in [TB] and applied 
to LEP data in m- As vacuum polarization in the 
space-like region is a smooth function of the squared 
momentum transfer, the accuracy of its determi¬ 
nation is only limited by the statistics and by the 
control of the systematics of the experiment. Also, 


1 





as at flavor factories the Bhabha cross section is 
strongly enhanced in the forward region, we will 
argue that a space-like determination of may 
not be limited by statistics and, although challeng¬ 
ing, may become competitive with standard results 
obtained with the dispersive approach via time-like 
data. 


2. Theoretical framework 


The leading-order hadronic contribution to the 
muon g-2 is given by the well-known formula Him 

a ds 

4 / - K{s) Imnhad(s + ^e), (1) 

^ do ^ 

where nhad(s) is the hadronic part of the photon 
vacuum polarization, e > 0, 


K{s) 



x^{l — x) 

+ {1 — x){s/m^) 


( 2 ) 


is a positive kernel function, and is the muon 
mass. As the total cross section for hadron produc¬ 
tion in low-energy e+e“ annihilations is related to 
the imaginary part of nhad(s) via the optical theo¬ 
rem, the dispersion integral in eq. Q is computed 
integrating experimental time-like (s > 0) data up 
to a certain value of s [amum. The high-energy 
tail of the integral is calculated using perturbative 
QCD [20] • 

Alternatively, if we exchange the x and s integra¬ 
tions in eq. Q we obtain m 

«r° = “/'^^(^-i)nhad[i(a;)], (3) 

^ Jo 

where nhad(i) = nhad(i) - nhad(O) and 

2 2 
X mf, 

t{x) = -^ < 0 (4) 

X — 1 


is a space-like squared four-momentum. If we invert 
eq. Q, we get x = {I — P) {t/2m'^), with /? = (1 — 
and from eq. (j^ we obtain 


a 


HLO 


a 

TT 


nhad(i) 



^ dt 

W 


(5) 


Equation © has been used for lattice QCD calcu¬ 
lations of [22]; while the results are not yet 
competitive with those obtained with the disper¬ 
sive approach via time-like data, their errors are 
expected to decrease significantly in the next few 
years [23j- 


The effective fine-structure constant at squared 
momentum transfer can be defined by 


a{q^) 


a 

1 — Aa(g^) ’ 


( 6 ) 


where Aa(g^) = —Ren(q^). The purely leptonic 
part, Aaiep((?^), can be calculated order-by-order 
in perturbation theory - it is known up to three 
loops in QED [24] (and up to four loops in specific 
q^ limits |2S]). As Imn(g^) = 0 for negative q^, 
eq. © can be rewritten in the form |26j 

= - / da; (1 - x) Aahad[t(a;)] • (7) 

Jo 

Equation ©> involving the hadronic contribution 
to the running of the effective fine-structure con¬ 
stant at space-like momenta, can be further formu¬ 
lated in terms of the Adler function [27], defined as 
the logarithmic derivative of the vacuum polariza¬ 
tion, which, in turn, can be calculated via a disper¬ 
sion relation with time-like hadroproduction data 
and perturbative QCD [23 [2H]- We will proceed 
differently, proposing to calculate eq. Q by mea¬ 
surements of the effective electromagnetic coupling 
in the space-like region (see also [9[)- 


3. AQ;had(^) from Bhabha scattering data 

The hadronic contribution to the running of a 
in the space-like region, Aohad {t ), can be extracted 
comparing Bhabha scattering data to Monte Carlo 
(MC) predictions. The LO Bhabha cross section 
receives contributions from t- and s-channel pho¬ 
ton exchange amplitudes. At NLO in QED, it is 
customary to distinguish corrections with an ad¬ 
ditional virtual photon or the emission of a real 
photon (photonic NLO) from those originated by 
the insertion of the vacuum polarization correc¬ 
tions into the LO photon propagator (VP). The 
latter goes formally beyond NLO when the Dyson 
resummed photon propagator is employed, which 
simply amounts to rescaling the a coupling in the 
LO s- and t-diagrams by the factor 1/(1 — Aa{q'^)) 
(see eq. Q). In MC codes, e.g. in BabaYaga [29] . 
VP corrections are also applied to photonic NLO 
diagrams, in order to account for a large part of 
the effect due to VP insertions in the NLO con¬ 
tributions. Beyond NLO accuracy, MC genera¬ 
tors consistently include also the exponentiation of 
(leading-log) QED corrections to provide a more 
realistic simulation of the process and to improve 
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the theoretical accuracy. We refer the reader to 
ref. [30] for an overview of the status of the most 
recent MC generators employed at flavor factories. 
We stress that, given the inclusive nature of the 
measurements, any contribution to vacuum polar¬ 
ization which is not explicitly subtracted by the 
MC generator will be part of the extracted Aa{q^). 
This could be the case, for example, of the contribu¬ 
tion of hadronic states including photons (which, al¬ 
though of higher order, are conventionally included 
in and that of W bosons or top quark pairs. 

Before entering the details of the extraction of 
Aahad(i) from Bhabha scattering data, let us con¬ 
sider a few simple points. In fig. (left) we plot the 
integrand (1—a;)AQ!had[i(2;)] ofeq. 0 using the out¬ 
put of the routine hadr5nl2 |3T] (which uses time¬ 
like hadroproduction data and perturbative QCD). 
The range x G (0,1) corresponds to t G (—oo,0), 
with a; = 0 for t = 0. The peak of the integrand 
occurs at Xpeak — 0.914 where fpeak — —0.108 
GeV^ and Aahad(ipeak) ^ 7.86 x 10“'* (see fig. 
(right)). Such relatively low t values can be ex¬ 
plored at e+e“ colliders with center-of-mass energy 
^/s around or below 10 GeV (the so called “flavor 
factories”) where 

i = ,8) 

6 is the electron scattering angle and rrie is the 
electron mass. Depending on s and 9, the inte¬ 
grand of eq. 0 can be measured in the range 
X G aJmax], as shown in fig. (left). Note 

that to span low x intervals, larger 9 ranges are 
needed as the collider energy decreases. In this 
respect, ^/s ^ 3 GeV appears to be very conve¬ 
nient, as an x interval [0.30,0.98] can be measured 
varying 9 between ^ 2° and 28°. It is also worth 
remarking that data collected at flavor factories, 
such as DA$NE (Frascati), VEPP-2000 (Novosi¬ 
birsk), BEPG-II (Beijing), PEP-II (SLAG) and Su- 
perKEKB (Tsukuba), and possibly at a future high- 
energy e+e“ collider, like FGG-ee (TLEP) [52| or 
ILG [33], can help to cover different and comple¬ 
mentary X regions. 

Furthermore, given the smoothness of the inte¬ 
grand, values outside the measured x interval may 
be interpolated with some theoretical input. In par¬ 
ticular, the region below Xmin will provide a rela¬ 
tively small contribution to while the region 

above Xm&x may be obtained by extrapolating the 
curve from Xmax to x = 1, where the integrand is 
null, or using perturbative QGD. 


The analytic dependence of the MG Bhabha pre¬ 
dictions on a{t) (and, in turn, on Aahad(i)) is not 
trivial, and a numerical procedure has to be devised 
to extract it from the dataQ In formulae, we have 
to find a function a{t) such that 


d(T 


dt 


data 


da 

dt 


(a{t),a{s)^ 


MC’ 


(9) 


where we explicitly kept apart the dependence on 
the time-like VP q;(s) because we are only inter¬ 
ested in a(t). We emphasise that, in our analy¬ 
sis, a{s) is an input parameter. Being the Bhabha 
cross section in the forward region dominated by 
the t-channel exchange diagram, we checked that 
the present q;(s) uncertainty induces in this region 
a relative error on the 9 distribution of less than 
^ 10“* (which is part of the systematic error). 

We propose to perform the numerical extraction 
of Aahad(i) from the Bhabha distribution of the t 
Mandelstam variable. The idea is to let a{t) vary 
in the MG sample around a reference value and 
choose, bin by bin in the t distribution, the value 
that minimizes the difference with data. The pro¬ 
cedure can be sketched as follows: 


1. choose a reference function returning the value 
of AQ;had(i) (and hence a{t)) to be used in the 
MG sample, we call it a{t); 

2. for each generated event, calculate N MG 
weights by rescaling a{t) —>■ a{t) + ^S{t), 
where i G [—iV, iV] and 6{t) is for example the 
error induced on a{t) by the error on Aahad(0- 
Being done on an event by event basis, the full 
dependence on a{t) of the MG differential cross 
section can be kept; 

3. for each bin j of the t distribution, compare 
the experimental differential cross section with 
the MG predictions and choose the ij which 
minimizes the difference; 

4. a{tj) -I- ^S{tj) will be the extracted value of 
a{tj) from data in the bin. Aahad(^j) can 
then be obtained through the relation between 
a{t) and Aahad)^)- 

We finally find, for each bin j of the t distribution. 


da 


dt 


j,da.ta 



^S{t),a{s)'j 


j,MC 


.( 10 ) 


^This was not the case for example in |16II17| : there a.{t) 
was extracted from Bhabha data in the very forward region 
at LEP, where the t channel diagrams are by far dominant 
and q(£) factorizes. 
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Figure 1: Left: The integrand (1 — x) Aa^ad[^(^)] 10^ a function of x and t. Right: AQ;had[^(^)] 10"^- 




e (deg) 9 {deg) 

Figure 2: Left: Ranges of x values as a function of the electron scattering angle 9 for three different center-of-mass energies. 
The horizontal line corresponds to x = a^peak — 0.914. Right: Bhabha differential cross section obtained with BabaYaga |29| as 
a function of 6 for the same three values of y/s in the angular range 2° < 9 < 90°. 
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We remark that the algorithm does not assume 
any simple dependence of the cross section on a{t), 
which can in fact be general, mixing s, t channels 
and higher order radiative corrections, relevant (or 
not) in different t domains. 

In order to test our procedure, we perform a 
pseudo-experiment: we generate pseudo-data using 
the parameterization of refs. [inilSl] and 

check if we can recover it by inserting in the MC 
the (independent) parameterization AQ!(,{jj(t) (cor¬ 
responding to a{t) of eq. 10 ) of ref. [3T] by means 
of the method described above. For this exercise, 
we use the generator BabaYaga in its most com¬ 
plete setup, generating events at = 1.02 GeV, 
requiring 10° < 9± < 170°, E± > 0.4 GeV and an 
acollinearity cut of 15°. We choose 6{t) to be the 
error induced on a{t) by the 1-cr error on Aahad(i)i 
which is returned by the routine of ref. [31] . we 
set N = 150, and we produce distributions with 
200 bins. We note that in the present exercise a(s) 
and all the radiative corrections both in the pseudo¬ 
data and in the MG samples are exactly the same, 
because we are interested in testing the algorithm 
rather than assessing the achievable accuracy, at 
least at this stage. 

In fig. 1^ is the result extracted with our 

algorithm, corresponding to the minimizing set of 
ij: the figure shows that our method is capable 
of recovering the underlying function AQ:had(i) in¬ 
serted into the “data”. As the difference between 
Aa^iad and Aa])^^ is hardly visible on an absolute 
scale, in fig. all the functions have been divided 
by Aa^i^j to display better the comparison between 
^“Ld and Aa°^([. 

In order to assess the achievable accuracy on 
Aahad(i) with the proposed method, we remark 
that the LO contribution to the cross section is 
quadratic in a(<), thus we have (see eq. §) 


1 Sa 6a 

77 — — 

2 (7 a 


( 11 ) 


Equation (11) relates the absolute error on Aohad 


with the relative error on the Bhabha cross section. 
From the theoretical point of view, the present ac¬ 
curacy of the MG predictions [30] is at the level 
of about 0.5%o, which implies that the precision 
that our method can, at best, set on Aahad(i) is 
(5Aahad(i) — 2 • 10““^. Any further improvement 
requires the inclusion of the NNLO QED correc¬ 
tions into the MC codes, which are at present not 
available (although not out of reach) [50] . 


From the experimental point of view, we remark 
that a measurement of from space-like data 
competitive with the current time-like evaluations 
would require an 0(1%) accuracy. Statistical con¬ 
siderations show that a 3% fractional accuracy on 
the integral can be obtained by sampling the 
integrand (1 — a:)Aahad[i(a:)] in ~ 10 points around 
the X peak with a fractional accuracy of 10%. Given 
the value of 0 ( 10 “^) for Aahad x = Xpeak, this 
implies that the cross section must be known with 
relative accuracy of ~ 2 x 10“^. Such a statistical 
accuracy, although challenging, can be obtained at 
flavor factories, as shown in fig. (right). With an 
integrated Inminosity of 0 ( 1 ), 0 ( 10 ), 0 ( 100 ) fb~^ 
at y6s = 1, 3 and 10 GeV, respectively, the angu¬ 
lar region of interest can be covered with a 0 . 01 % 
accuracy per degree. The experimental systematic 
error must match the same level of accuracy. 

A fraction of the experimental systematic error 
comes from the knowledge of the machine luminos¬ 
ity, which is normalized by calculating a theoretical 
cross section in principle not depending on Aohad- 
We devise two possible options for the normaliza¬ 
tion process: 

1 . using the e^e“ —>■ 77 process, which has no 
dependence on Aohad, at least up to NNLO 
order; 

2. using the Bhabha process at t ~ 10“^ GeV^ 
{x ^ 0.3), where the dependence on Aahad is 
of 0 ( 10 “®) and can be safely neglected. 

Both processes have advantages and disadvantages; 
a dedicated study of the optimal choice goes beyond 
the scope of this paper and will be considered in a 
future detailed analysis. 


4. Conclusions 


We presented a novel approach to determine the 
leading hadronic correction to the muon g-2 using 
measurements of the running of a{t) in the space¬ 
like region from Bhabha scattering data. Although 
challenging, we argued that this alternative deter¬ 
mination may become feasible using data collected 
at present flavor factories and possibly also at a fu¬ 
ture high-energy e+e” collider. The proposed de¬ 
termination can become competitive with the ac¬ 
curacy of the present results obtained with the dis¬ 
persive approach via time-like data. 
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Figure 3: The extracted function compared to the function used in the pseudo-data (see text). The 

functions Aa^^^(t) it S(t) are shown to display the range spanned by the MC samples. All functions have been divided by 
Aa^^d(i)- The tiny difference between Aa^ad binning discretization. 
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